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WITH EXPERIMENTAL DATA 

By John N. B. Livingood and Anthony J. Diaguila 


SUMMARY 

A procedure is presented for the rapid determination of the core 
dimensions of a heat exchanger having one dominating film resistance. The 
length of the exchanger in the direction of the primary fluid flow and 
the Reynolds number of this flow are determined graphically from three 
trail solutions of the heat-flow and pressure-drop equations. Methods for 
determining the other two dimensions are also discussed. 

By use of experimental data obtained for a shell and tube liquid- 
metal-to-air heat exchanger, the calculation procedure presented herein 
is verified. Results are within the accuracy of the spread of the exper- 
imental data. The use of approximate flow conditions yields adequate core 
dimensions in the examples given herein. 


INTRODUCTION 

The use of heat exchangers in high-speed, high-altitude aircraft is 
receiving more and more attention. For such applications, the size and 
weight of the heat exchanger and the power required to drive the coolants 
through the exchanger become the predominating factors . In view of this 
fact, a method for optimizing any one of the heat -exchanger parameters 
(power, weight, volume, or frontal area) against any other one was devel- 
oped and reported in reference 1. This optimization was determined for 
heat exchangers with one dominating film resistance and included seven 
typical configurations of reference 2. 

Recently, an attempt was made to investigate the feasibility of a 
gas-to-gas heat exchanger for use in reducing the temperature of compres- 
sor bleed air prior to its use as the turbine coolant in high-speed, 
high-altitude turbojet engines. For such an application, the inlet and 
exit states of both fluids, the available pressure drops of both fluids, 
the temperature change of one fluid, and the heat capacities of both fluids 
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may be prescribed. For prescribed pressure drops of both fluids, a con- 
ventional calculation procedure for determining the heat-exchanger size 
becomes lengthy and involved. A method which determines a gas-to-gas 
heat -exchanger core size with minimum time and effort and which is suit- 
able for the aforementioned application was developed and presented in 
reference 3 for a prescribed core configuration and a single set of fluid 
conditions. This procedure has since been generalized for any gas-to-gas 
heat -exchanger core configuration and a range of fluid conditions in 
reference 4. The sizes and weights of a number of heat -exchanger cores 
for possible use in aircraft flying at Mach 2.5 and 70,000 feet were 
determined by the method of reference 3 and presented in reference 5. 

In nuclear reactors, and in some possible aircraft engines, heat 
exchangers that employ a liquid or liquid metal are of importance. Such 
exchangers have one dominating film resistance and may be optimized by 
the method of reference 1. For this type of heat exchanger, the heat ex- 
changed depends essentially on the conditions of the gas (primary coolant). 
Consequently, only the pressure drop of the primary fluid is considered 
in this application. The core dimensions of this type of heat exchanger 
can be determined by use of a modification of the method presented in 
reference 3. This report presents this modified method and compares, for 
specified conditions, dimensions determined by use of this method with 
those of an experimental shell and tube liquid-metal-to-air heat exchanger 
(ref. 6). The investigation was made at the NACA Lewis laboratory. 


SYMBOLS 

The following symbols with consistent units are used: 

I 

A heat -transfer area 

A' free-flow area 

Ap frontal area 

Cp specific heat at constant pressure 
D inside diameter of heat -exchanger shell 

d hydraulic diameter 

f friction factor 

g acceleration due to gravity 

h heat-transfer coefficient 



NACA RM E56K26a 


3 


k thermal conductivity 

L heat -exchanger core length 

l half fin length 

Pr Prandtl number 

p pressure 

R gas constant 

Re Reynolds number, wd/A'p 

St Stanton number, hA'/wcp 

T temperature 

Tu heat-transfer parameter (number of transfer units, denoted as NTU 
in ref. 2) 

tp fin thickness 

U over-all heat-transfer coefficient 

V velocity 

v specific volume 

w weight -flow rate 

a heat -transfer surface area per unit volume 

Tip fin effectiveness 

tyj, thermal effectiveness (denoted as 6 in ref. 2) 
tJq surface effectiveness 

viscosity based on film temperature 
p density 

a ratio of free-flow to frontal area, A'./Ap 

Subscripts: 

ex exit 
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f fin 

i inlet 

max maximum 

min minimum 

n no-flow direction 

1 heat exchanger side and fluid with finite heat resistance 

2 heat exchanger side and fluid with negligible heat resistance 


DEVELOPMENT OF METHOD FOR CALCULATING CORE DIMENSIONS 
Heat -Exchanger Equations 

Under the assumption that one heat resistance is negligible, that 
is U = T] 0 h 1 , the number of transfer units in the heat exchanger may be 
written as (ref. l) 


HphjAi 

" W 1 C P>1 


( 1 ) 


When the heat transfer coefficient h L is replaced by its .equivalent 
(eq. (14) of ref. 3), 


h l * (ReSt)! 


( 2 ) 


the area A 1 is replaced by its equivalent (eq. (15) of ref. 3), 


a l 

A 1 - ^ A i L l 


( 3 ) 


and w-l is expressed in terms of the Reynolds number, equation (l) may 
be written 


Re^Tuc! 

Ll “ ^ItRistTI 


( 4 ) 


where L x is the exchanger core length on the fluid side with the finite 
heat resistance. 
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A second expression for the same length Lp is obtainable from the 
pressure-drop equation (eq. (17) of ref. 3 with end losses neglected), 
which may be written 


[ 2 g 

gd-1 API 1 + PI v l,ex ~ v l,i 

Re l li l( v ljex + v l,i) ’ 2 + Vl - 1 . 


( 5 ) 


where 


RTi i 

v-i ,• = *— 

1,1 Pl,i 


_ R(Tl,i + gTp) 
V1 ' ex Pl,i - 4?1 


and p. is evaluated at the film temperature. 


For prescribed values of pp^, Tp^p, ^1* w l c p,l^ w 2 c p 2> an< ^ 

Tg^p , and for a prescribed core configuration, equations (4) and (5) 
become a pair of equations in two unknowns. Rep and Lp. A method of 

solution for these equations will be given later. For the conditions 
previously prescribed, values of the other dimensions Lg and 1^ of 

the heat -exchanger core can be determined in the following way. With the 
values of Rep and h. obtained from the solution of equations (4) and 

and (5), the frontal area for the primary coolant side (or the product 
LgL n ) can be determined from the continuity equation for the primary fluid, 

that is, 


wp = apRep L 2 Ln (6) 

If either Lg or 1^ is known or calculable, the other length (i^ or 
Lg) can be determined from equation (6) . Three cases will be considered. 

Case 1. - For somek applications of the type of heat exchanger con- 
sidered herein, installation or other considerations may require that 
either L2 or Ln be restricted in size. In this case, the other length 
can be determined directly from equation (6) . 

Case 2 . - In other applications, it is conceivable that the velocity 
of the secondary fluid may be restricted to a certain value. In this 
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case, the value of Lq can be determined from the continuity equation 
of the secondary fluid written in the following form: 

w 2 “ P2 v 2° r 2 L l L n ( 7 ) 

With the appropriate substitutions, L2 can then be found from equa- 
tion (6) . 

Case 5« - If neither of the limitations of cases 1 or 2 apply, it 
may be necessary to assume a length for either L 2 or L n , and then 

obtain a length for the other dimension (i^ or L 2 ) from equation (6). 

In this way a series of heat -exchanger geometries can be obtained, and 
the particular selection is left to the designer. 


Calculation Procedure for Solving Equations (4) and (5) 

From the prescribed conditions stated previously and from the heat- 
flow equation 


*1^,1 ^1 + w 2°p,2 ^2 - 0 


( 8 ) 


the value of ZJ£T 2 is obtained, 
temperature . The equation 


The values . of c~. are based on bulk 


m. 


max 


T l,i - T 2,i 


0 ) 


may then be solved for rjip. For various flow conditions, reference 2 
presents plots of against Tu with (vc p ) min /(wc p ) max (or 

w l c p,l/ w 2 c p,2 ^ or cases considered herein) as parameter. From the 

prescribed conditions, equation (9), and the appropriate curve in refer- 
ence 2, the value of Tu is then determined. 


The values of and cr^ are obtainable from the prescribed core 

configuration) and t) Q for either side is obtained from the following 
expressions (eqs. (12) and (13) of ref. 3): 


no 


Af» 

1 - x (1 " nf) 


( 10 ) 
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where 

v-4P****^» (-) 

Friction and heat -transfer data [fp and (ReSt)p for use in equations 

(5) and (4)] in terms of the Reynolds numbers are obtainable from refer- 
ence 2 for various core configurations or from equations applicable to 
the type of passages considered in the selected core. 

From the preceding information, equations (4) and (5) can each be 
solved for Lp for a series of assumed values of Rep. The intersection 

of the curves representing the corresponding values of Rep and Lp from 
the two equations yields the desired correct values of Rep and Lp. 

VERIFICATION OF METHOD BY USE OF EXPERIMENTALLY DETERMINED VALUES OF Tu 

The method for calculating the core dimensions of a heat exchanger 
presented herein may be verified with the aid of experimental data ob- 
tained for a liquid-metal-to-air shell and tube, heat exchanger (ref. 6). 

A schematic diagram of this experimental exchanger is shown in figure 1. 
Air flowed through 241 tubes of 3/l6-inch outer diameter, 0.016-inch wall, 
and 28-inch length (L/d = 180.). A 4.25-inch-inside-diameter shell sur- 
rounded the tube bundle. Sodium flowed over the tubes as indicated in 

figure 1. For this exchanger, ap is 99.26 feet -1 and ap is 0.3205. 
Experimental data necessary for the verification are presented in table I. 

Friction and heat -transfer data were correlated in reference 6 by 
use of appropriate and well-established correlations . These same cor- 
relations will be used herein. The conventional single -tube heat- transfer 
correlation corrected for an L/d ratio of 180 and with an assumed value 

of Pr 2 / 3 equal to 0.75 will be employed; this correlation is 

(ReSt ) p = 0.028 Rep’ 8 (12) 

The Karman-Nikuradse friction correlation 



is also used. 
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For any set of experimental data, the value of (ReSt)^ can be ob- 
tained from equation (.12) for the experimentally determined Re^. With 
these values, Re^ and (ReSt)^, the experimental exchanger length 
of 28 inches, the values of a-^ and cx^ for this exchanger, and tjq 

of 1 (there are no fins), the correct value of Tu for this set of data 
can be obtained from equation (4) . The verification of the method of 
calculating the heat-exchanger core dimensions presented herein can be 
demonstrated by use of these experimentally determined values of Tu as 
follows . 

For any set of data, and the corresponding experimental value of 
Tu, equation (4) can be solved for for three assumed values of Re-^. 

The values of (ReSt)^ corresponding to the assumed values of Re^ are 
obtained from equation (12) . For the same three assumed values of Rej_, 
values of f^ are obtained from equation (.13) . With these values and 
the experimentally determined values of p-j^, Ap-|_, and AT^, three 

solutions of equation (5) can be determined. The intersection of the 
curves representing the solutions of equations (4) and (5) yields the 
desired values of and Re^. Any discrepancy between this value of 

L 1 and the design length of 28 inches must therefore be attributed to 

(l) the discrepancies between the heat-transfer and friction equations 
(eqs. (10) and (ll)) and the measured heat-transfer and friction data, 
and (2) the differences in the heat flow between the two fluids (about 
8 %, according to ref. 6)-. 

Since the experimental exchanger is of the tube-shell type, the 
frontal area on the primary fluid side is circular in shape, and equation 
(6) may therefore be written 

wi = a 1 Re 1 ^ (14) 

where the frontal area Lgl^ is replaced by «D 2 /4. Equation (14) can 

be solved for D, and these values of D can be compared with the design 
dimension of 4.25 inches. 


APPLICATION OF METHOD FOR DESIGN PURPOSES (WHEN 

EXPERIMENTAL DATA ARE UNAVAILABLE) 

When one is faced with the problem of designing a heat exchanger 
and experimentally determined values of Tu are not available, it is 
necessary to obtain values of Tu by use of equation (9) and an 
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appropriate relation between ( wc p)min/ ( wc p)max ( e< J ual to 

w l c p i/ w 2 c p 2 herein), and Tu for the particular flow conditions in- 
volved. If such flow conditions are well defined (such as crossflow, 
counterflow, or parallel flow) , accurate values of Tu are available from 
reference 2. If flow conditions are not well defined, it may be necessary 
to approximate the values of Tu by use of certain assumed flow condi- 
tions. For the experimental heat exchanger previously described, crossflow- 
counterflow conditions prevail. Since reference 2 does not contain a 
Tu against T]ip plot for this flow condition, it was necessary to assume 
an approximating condition; crossflow was assumed. When design conditions 
identical to the inlet and exit conditions of the experimental exchanger 
are considered and when values of Tu for assumed crossflow are applied, 
approximate core dimensions can be determined. These approximate dimen- 
sions are compared with those of the experimental exchanger. It should 
be emphasized that the calculation procedure is applicable for any 
chosen set of inlet and exit conditions . The use of experimentally meas- 
used conditions herein is made solely for the purpose of comparison. In 
this way, it is possible to determine whether the use of approximate flow 
conditions results in calculated dimensions close to those obtained with 
the use of true flow conditions. 

From the experimental inlet and exit conditions and the assumed 
crossflow condition, values of Tu for each set of data can be obtained. 
Values of qip are found from equation (9). For these typ values and 

the corresponding values of (wcp) min /(wcp) max (determined from the infor- 
mation presented in table i), the Tu values are read from figure 2; 
figure 2 is reproduced from reference 2 and applies for crossflow condi- 
tions. For three assumed values of Re-j_, and the corresponding values of 

(ReSt)i and f^ obtained from equations (12) and (.13), three solutions 
for Lp of each of equations (4) and (5) can be determined for each set 

of data in table I. The intersection of the curves through the three 
solutions of each of equations (4) and (5) gives the desired values of 
and Re^ . 


RESULTS AND DISCUSSION 

Initial calculations were made according to the procedure discussed 
in the section VERIFICATION OF METHOD BY USE OF EXPERIMENTALLY DETERMINED 
VALUES OF Tu. The results of these calculations are as follows: 
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a 

a 

Pd 

Experi- 
mentally 
deter- 
mined Tu 

Re ^ 

Ll, 

in. 

in. 

*D 2 /4, 
sq in. 

Volume, 
cu in. 

1 

2.58 

31,750 

28.6 

4.27 

14.3 

408 

2 

3.05 

14,950 

29.1 

4.11 

13.3 

387 

3 

2.72 

25,700 

28.8 

4.24 

14.1 

406 

. 4 

2.48 

33,000 

27.8’ 

4.51 

16.0 

445 

5 

2.67 

24,700 

28.2 

4.37 

15.0 

424 

6 

2.71 

22,000 

27.7 

4.36 

14.9 

413 

7 

3.07 

13,400 

28.6 

4.10 

13.2 

377 

8 

2.33 

43,000 

27.5 

4.54 

16.2 

446 

9 

2.29 

45,300 

27.4 

4.58 

16.5 

451 


Figure 3(a) shows the solutions of equations (4) and (5) for the 
assumed values of Re^, and the intersection point of the curves joining 
these solutions (L]_ and Re^.) for run 5. Comparison of the average 
calculated values of L^, nD^/4, and volume for the nine runs (28.2 in., 

14.8 sq in./ and 417 cu in., respectively) with those of the experimental 
exchanger (28 in., 14.2 sq in., and 397 cu in., respectively) resulted 
in discrepancies of about 1, 4, and 5 percent, respectively. For some 
runs, these discrepancies increased to as much as about 4, 16, and 13 
percent, respectively. Reference 6 shows that the measured data deviated 
from the heat-transfer and friction correlation equations by as much as 
15 percent. Since the correlation equations were used in the calculations 
just discussed, it may be concluded that the calculation method presented 
herein is verified by giving results within the accuracy of the experimen- 
tal data . 

Calculations were also made by the procedure described in APPLICATION 
OF METHOD FOR DESIGN PURPOSES (WHEN EXPERIMENTAL DATA ARE UNAVAILABLE) . 

As stated previously, inlet and exit conditions identical to the experi- 
mental values listed in table I were selected, but crossflow was assumed 
as an approximation for the flow conditions within the exchanger. The 
calculations, which for this case may be termed approximate because of 
the assumed crossflow, were made for the same trial values of Re-^ that 

were assumed for the other calculations . The results are 
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Run 

Cross- 

flow, 

Tu 

Re-^ 

*L> 

in. 

D, 

in. 

*D 2 /4, 
sq in. 

Volume, 
cu in. 

1 

2.40 

31,900 

28.1 

4.26 

14.1 

400 

2 

2.52 

16,100 

24.7 

3.96 

12.3 

305 

3 

2.30 

27,800 

24.9 

4.02 

12.7 

316 

4 

2.18 

34,400 

24.8 

4.41 

15.3 

381 

5 

2.41 

25,600 

25.7 

4.30 

14.5 

372 

6 

2.65 

22,200 

27.5 

4.34 

14.7 

406 

7 

2.85 

13,800 

26.9 

4.03 

12.8 

345 

8 

2.16 

44,200 

25.6 

4.49 

15.8 

403 

9 

2.12 

46,500 

25.6 

4.52 

16.0 

411 . 


Figure 3(b) presents the solutions of equations (4) and (5) for the 
assumed values of Re^ and the intersection point of the curves Joining 

these solutions for run 5. Average values of L^, jtD^/ 4, and volume for 

the nine runs (26 in., 14.2 sq in., and 371 cu in., respectively) now differ 
from the design values by about 7, 0, and 6 percent, respectively. For some 
runs, discrepancies increased to as much as 12, 13, and 23 percent. 

From the calculations presented, it may be concluded that the use of 
approximate flow conditions yields adequate core dimensions. For precise 
calculations, an accurate knowledge of flow conditions is essential. 


SUMMARY OF RESULTS 

The results of this investigation are summarized as follows: 

1. A method is presented for the rapid determination of the core 
dimensions of a heat exchanger having one dominating film resistance. 

Three trial solutions of the heat-flow and pressure-drop equations are 
sufficient for determining the heat -exchanger length in the direction of 
the primary fluid flow and the Reynolds number of this flow. Methods for 
determining the other two core dimensions are also discussed. 

2. The method is verified with experimental results obtained from a 
shell and tube liquid-metal-to-air heat exchanger. For experimentally 
determined values of the heat -transfer parameter, the average values of 
the exchanger core length (on the fluid side with the finite heat re- 
sistance), frontal area, and volume differed from the experimental ex- 
changer values by about 1, 4, and 5 percent, respectively. When the heat- 
transfer parameter values were found from an available heat-transfer 

chart for flow conditions that approximated those existing in the exchanger, 
the average values of the exchanger core length (on the fluid side with 
the finite heat resistance), frontal area, and volume differed from the 
experimental exchanger values by about 7, 0, and 6 percent, respectively. 
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3. The discrepancies between the calculated dimensions and those of 
the experimental heat exchanger, when utilizing the experimentally deter- 
mined values of the heat-transfer parameter, result from the scatter in 
the heat-transfer and friction data as well as the apparent differences 
in the heat flow between the two fluids . 

4. For the calculations considered herein, the use of approximate 
flow conditions gave adequate core dimensions. An accurate knowledge of 
flow conditions is essential for the determination of precise core 
dimensions . 


Lewis Flight Propulsion Laboratory 

National Advisory Committee for Aeronautics 
Cleveland, Ohio, November 23, 1956 
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Figure 2. - Performance of crossflow heat exchanger with fluids 
unmixed (ref. 2). 











